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One of the simplified solutions of the Dirac equations with the pure Coulomb potential 
given in a paper published in 1985 is pseudo. The original paper solved the Dirac equa- 
tions by introducing a transformation of functions with two strange parameters a and h to 
transform the original system of the first-order differential equations into two uncoupled dif- 
ferential equations of second order. However, not only the given eigenvalues sets violate the 
uniqueness of solution but also the said second-order equations are not any necessarily math- 
ematical deduction. In order to determine the introduced parameters, formally, the author 
actually introduced some self-contradictory mathematical formulas, such as sinh 9 = 2ab, 
cosh 61 = + b'^, tanhO = -Za/k, a"^ - ^ I, b ^ sinh(6'/2) and a = cosh(6'/2). But one 
has not known the value of the parameters a and b all the while, whereas the parameters 
were insensibly deleted in the given second-order Dirac-Coulomb equation last. One can- 
not recover any result given in the paper by making corresponding correctly mathematical 
calculations. 
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I. INTRODUCTION 

It is well known that the exact solutions of the wave equations[ll] are very important. Because 
the judgment on whether or not the approximately solution of the wave equations exist usually 
dependents on the existence of exact solutions and some real laws are often obtained from 
the exact solution of the wave equations. For the Dirac equation for a single particle in the 
Coulomb field, quantum mechanics textbooks generally adopt the exact solution that was first 
given by Darwin and Gordon[^Q]. Of course, there are some mathematical problems needing 
to discuss in the Darwin-Gordon solutionjsl. However, some other exact solutions of the second- 
order Dirac-Coulomb equation appear simphfied but are clearly incorrect in characterflfl. In 
form, it derived those so-called simplified solution by first transforming the system of the first- 
order Driac-Coulomb equation into the Schrodinger-like or the Klein-Gordon-like equations then 
solving the corresponding second-order differential equations to write the distinguished energy 
eigenvalues. Nevertheless, one cannot recover the corresponding solution basing on the introduced 
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mathematical methods in the original papers. It has been pointed out that many given solutions 
and the corresponding formulas of the energy levels are actually not the necessary mathematical 
deductions of the said second-order Dirac-Coulomb equation. In addition, it is worse that some 
so-called second-order Dirac-Coulomb equations given in those published papers can not be yielded 
from the original system of the Dirac- equations of first-order in the Coulomb field. They are not 
any necessary mathematical deduction of physics and mathematics yet. 

Here we show that a historically formal simplified solution of the Dirac-Coulomb equation given 
in a paperQ] published in Physical Review 22 years ago is a pseudo solution. For the said second- 
order Dirac-Coulomb equations in which two equations were written in the same form by using 
sign "lb", two eigenvalues set should be given and they are actually different. It is well known that 
two different sets of the energy eigenvalues for the same quantum system violate the uniqueness of 
solution. However, in the original paper, only one of the eigenvalues sets was given and the other 
was thrown out of all reason. This case still exists in other papers today [? ][? In particular, 

in order to write the second-order Dirac-Coulomb equation, the author introduced two strange 
parameters a and b and afterward gave some self-contradictory expressions such as sinhS = 2ab, 
cosh^ = a'^ + b'^, tanh 9 = —Za/k, — b'^ = 1, b = sinh (6/2) and a = cosh (9/2). However, one has 
not known the values of parametera and b. By using the mathematical method introduced therein 
or by using other correct mathematical methods, one cannot recover any so-called Schrodinger- 
like or Klein- Gordon-like equations given by the author. Consequently, in the mentioned paper, 
the claimed simplified solutions of the Dirac-Coulomb equation, the corresponding formula for the 
energy levels and the corresponding mathematical procedures are pseudo. 

II. ORIGINAL FORMAL SIMPLIFIED SOLUTION OF DIRAC-COULOMB EQUATION 

Many authors claimed that they obtain simplified solution of the Driac-Coulomb equation. In 
1985, Su considered that, with the use of a simple similarity transformation which brought the 
radial wave equations of the Dirac-Coulomb problem into a form nearly identical to those of the 
Schrodinger and Klein-Gordon equations, he derived simplified solutions to the Dirac-Coulomb 
equation for both the bound and continuum states following the familiar standard procedurfe 
adopted in the derivation of the conventional solutions. He considered that to obtain the desired 
form of the second-order radial equations he could still work with a first-order partial differential 
equation rather than with the second-order Dirac equation widely employed in the derivation of the 
simplified solutions, and thus he can avoid the task of reducing the solutions of the second-order 
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equations to those of the original Dirac equations. Here we only check the mathematical procedure 
for deducing the eneregy eigensolutions for the bound state in the original paper. The author first 
wrote the radial Dirac-Coulomb equation in the following form 

(1) 





with 
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(2) 



where 



(3) 



mc2 cosh^ + he [sinhe (|r + ^) - ^] 
:- {mc2 sinh^ + ;ic [cosh (I: + i) - f ] } 
mc^ sinh 6 + he [cosh ^' + ^) + 7] 
: - {mc2 cosh^ + he [sinh^ + f ) + f^] } 
and A; = lb + i) , a = 1^ being the fine-structure constant. The strange parameters a and h were 
claimed the real constants by the author. Then it was introduced that 



A'-- 
B' 
C" 
D' 



a^ + b'' . 2ab 

COSh6'=-:5 rrj, Smh 6/ = -jj -Tj. 

— 0^ — 0^ 
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and it was selected that 

a — 6 = 1, a = cosh I 2 ) ' ^ ~ ^^'^ 
with At the same time, the author also introduced the other function 

he = -Za/K 

It was alleged that one attained great simplification in solving the radial equations 
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where = ^1, 7 = (i + 5) - Z^oi^ 
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Finally, without some mathematical calculations, it was 



also alleged that one should derive the equation 
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E'^ - m^c^ ^ 2EZa\ 7^ ± 0)7 
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and the distinguished Dirac formula of the energy levels in the Coulomb field for bound state 
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III. ORIGINAL EIGENSOLUTIONS SETS VIOLATE THE UNIQUENESS OF 

SOLUTION 

In fact, the system of the second-order equations is not always equivalent to the corresponding 
first-order differential equations. Firstly, we don not know what it means for the parameters a and 
6and how to eliminate the two parameters to derive the second-order equations ([8|) from those new 
definitions such as from ([I]) to ([7]). We notice that the equation ([8]) given in the original article 
should include four equations as follows 



l_d / 2^) + / E^^rp^ + - ) R = 

dr \ dr J \ h/c^ her / 

(10) 



1 A /'^2dfi\ I t E^~m^d^ , 2EZa _ 7_±7 \ R _ n 
T^dr V dr ) ^ \ I" her | n — U 



1 _d (^2dQ\ , ( E^-m'^c'^ , 2EZa _ 7-7 \ D — (\ 
dr y dr J ^ \ '~ her j W — 'J 

1 _d_ (^2dQ\ _|_ ( E^-rn'^e'^ , 2EZo 
r^ dr \ dr J \ hFc? her 

Each of the equations has its own eigenvalus and eigensolutions set, and these eigenvalues are 
usually different from each other. It is incorrect for giving only one of the formulas of the energy 
levels and alleging to recover the distinguished Dirac formula in mathematical and physical logic. 
One easily finds the eigenvalue sect of the general equations 

1 d f 2dA\ , [E'^-m'^c^ ^ 2EZa 7^ =F 7 



\d (2dQ\ ( EWcl + 2|Za _ 7!+7 ) Q = 
■ dr \ dr I ' \ hf^r-^ ' nr.r / ^ 



r"^ dr\ dr ) ^ \ K^c^ ~'~ her J ^ ^ 

Because some details are uselly ignored by using the corresponding special function, we directly 
solve every differential equation to find the energy eigenvalues. By introducing the substitution 

(12) 

The equation ([TT|) becomes 

d^M ( m'^c^-E'^ ^Za^ 72 

+ ^27^ + 1— ^-^^W = (13) 



dr^ y ItP'c'^ her 
This equation has the asymptotic solutions with E < mc^ satisfying the boundary condition 




M = e^p\-\l r (r ^ 00) (14) 



It is assumed that the formal solution of the equation (jl3p takes the form 




M = exp(-^/ ^^^^ r\u (15) 
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We have 



dM _ 
dr — ^^P 



r -T- 
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d^u 
dr^ 



m^c'^—E^ du _j_ m^c' 



(16) 



Substituting into the equation (jl3p . we obtain 



(i^u ^ m'^c^ - du ^ /2EZa 7^T7\ _q 



Now, seeking the power series solution, let 



oo 
n=0 



s+n 



(17) 
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Hence 



^ = f; (s + n) d„r^+"-i, ^ = ^ (s + n) (s + n - 1) d„r^+"- 



(19) 



n=0 n=0 

Substitute and (fTUj) into the equations PT]) . we have 



n=0 



[(s + n) (s + n-l)-(72^7)]4 
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It gives the recursive relation of the coefficients of the power series 
[(s + n) (s + n- 1) - (72^7)] fi„ 



2a/ TTT^^ (s + n-1) 



he 



dn-l = (21) 



The power series (llSp naturally give the initial value condition: = d_2 = • • • = and do 7^ 0. 
Putting let n = and substituting it into the recursive relation ()2ip reads s (s — 1) — (7^ + 7) = 0, 
it gives 

1 ± (27 + 1) 



(22) 



Thus we have the multi- values of s for the four equations in pUj) respectively 



7 

1-7 



7 

1-7 



■5Q2 



1+7 

-7 



(23) 



1 + 7 

-7 

Since the wave function has to be normalizable we must choose the value of s to be more thanl 
but not only positive sign. In form, for the solutions ([23]), we can but choose 



SRi = 7> SR^ = l+ 7, SQ, =7, SQ2 = 1 + 7 



(24) 
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Combining the expression p2|) . p5|) . (jT8|) . (p4|) . and making using of the value of 7, the com- 
pletely formal solution of the equations (jlOp can be written as follows 

R, = exp (.y^^g^r) gd„rV0-+i/2)^-^^"^+n-i 

m^d^-E-i ^ g (i„rV(i+l/2)'-^2„2+„ 





= exp ^- 




= exp 


Qi 


= exp ^- 


Q2 


= exp ( - 



"=o (25) 



—A g d„,rV(i+l/2)"-22a2+n-l 
/ n=0 

/ n=0 



When j = 1/2, 7 = Vl — ^^a^ < 1, the first and the third expression are divergent at the origin 
of the coordinate system 



lim^i?i = lim^exp ( -^/™?g^r] g d„rV(i+i/2)^-z2a2+n-i ^ ^ 



} ' <;"=o ^ (26) 

limQi = lim exp ( I V d„rV(i+i/2)'~z2a2+n-i ^ ^ 

r->0 r^O \^ V c y 

implying that the first equation and the third equation in (llOp have no solution which satisfy the 
boundary conditions. We know that the above divergence have been called "mild divergence" by 
someone. In (j25p . the second and the forth expression are finite at the origin of the coordinate 
system 



hm R2 = lim exp (-.[^^^A g d„rV(J+i/2)'-^^"2+n = q 

r^O r^O V ^ / n=0 

hm Q2 = hm exp (-\ h^t~i'^\ \ g d„rV(i+i/2)'-2^«^+" ^ g 



(27) 



implying that the second and the forth equation in (|10p seem to have the significative solutions. 
Comparing (j26p and (j27p constructs the first kind of contradiction to the second-order differential 
equation ([TO]) . 

On the other hand, if we accept the subjective definition of the so-called "mild divergence" 
or "weak divergence" we would have two eigenvalues set corresponding to the formula of energy 
levels in the Coulomb field. Form (j25p . combining the first and the third expression in one form 
denoted by Ai and combining the second and the forth expression in another form denoted by A2 
respectively yield 

Ai = exp (-M^A g d„rV0-+i/2)^-^^a2+„-i 

> , < "=o ^ (28) 

A2 = exp i-Jrri^r] E d„rVa+i/2)^-^^«Hn 



\ ■ / n=0 

In order that the wave functions remain normalizable we must require that the series for u so the 
recursive relation ()2ip terminate at any term with the power rir, that means that dn^ 7^ and 
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drir+i = = • • • = 0. Substituting n = n^+i into ([2T]) . we have 



2V^^(. + n.)-^ = (29) 



So that we obtain the Dirac formula in form 



{rir+s) 

According to (I24p , we finally obtain the eigenvalues of the energy 

^2 



£■ = , (30) 

-L + — TT 



1+- 



^ ^ , ^ (31) 



^ (^„,, + l + -y/(j + l/2)2-z2„2^ 

It implies that the solutions of the equations ^ given in the original paper violate the uniqueness 
of solution of the wave equations. This is the serious mathematical and physical mistakes. In 
addition, the formula ([9]) is not a necessarily mathematical and physical deduction. It constructs 
the second kind of contradiction to the second-order differential equation (jlOp . 



IV. THE SECOND-ORDER DIRAC EQUATION ^ IS PSEUDO 

In fact, one cannot transform the original radial Dirac-Coulomb equations with first-order into 
the Schrodinger-like equation or the Klein- Gordon-like equation with the second-order ([8]). It means 
that the equation is not equivalent to the original system of first-order Dirac-Coulomb equations. 
Trying to explain the meaning of the definition such as sinh^, cosh0 and tanh^ make us puzzled, 
the expressions dH, ([5]) and ([6]) read 

sinh^ = 2ab, cosh0 = a'^ + 6^, tanh^ = —Za/k 

' (32) 

0^-62 = 1^ a = cosh (0/2) , 6 = sinh(6l/2) 
These expressions given in the original paper early or late are in contradiction with each other. 
Why did not the author direct give the value of a and 6, but deleted them finally expression like 
magic without any mathematical operation?. Now writing the equations ([T]) with ([2]) and ([3]) in 
the separate form as follows 

{mc2 cosh B^Uc [sinh 6' (|: + ^) - 4^] } ^ 

- {mc2 sinh B^hc [cosh (|- + 1) - ^] } Q = ^^^^ 
{mc2 sinh e + hc [cosh (|, + 1) + ^] } i? 

- {mc2 cosh B + Uc [sinh (-| + 1) + f^] } Q = 
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By using the definition it becomes 



2abhc dR , 
a^-b'^ dr 



(a2+62) 
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2abhc dQ _ 
a'^—b'^ dr 



+ 
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(34) 



Q = 



Because it is not known for value of the parameter a and 6, we cannot obtain the formula of the 



energy levels in Coulomb field. According to the conflicting expressions 
it can only yield 



by using a? — iP' = 1, 



(35) 



2ahhc^ + [{a^ + 6^) mc^ - E + {2ab - Za) ^] R 
- («' + - [2a^rnc2 + {a^ + - k) ^] Q = 

(a2 + 62) ncf + [2abmc^ + {a^ + b^ + k)!f]R 
-2abhc^ -[{a^ + 6^) mc^ + E + {2ab -Za)^]Q = 

This is the original shape of the introduced expression dT]) . It is clear that this system of differential 
equations are not equivalent to the original Dirac-Coulomb equations, and one cannot translate 
them into the so-called Schrodinger-like equations ([8]). Only when one knows the value of the 
parameters a and b can translate the system of first-order differential into the corresponding second- 
order differential equations without the undetermined parameters. Consequently, the second-order 
equations ([8]) are the pseudo equations, and the formula of the energy levels ([9]) is not a necessary 
mathematical deduction. 



V. CONCLUSIONS 



In the present paper, we have shown that the original paper published in 1985 used many 
self-contradictor definition expressions to finally written the so-called simplified solutions of the 
Dirac-Coulomb equation and all of the given results in the original paper cannot be recovered via 
the strict mathematical calculating. Such kind of simplified solutions is the pseudo solution of 
the Dirac equation in the Coulomb field. The corresponding second-order Dirac-Coulomb equa- 
tions are the pseudo second-order Dirac equation. In fact, it is very simple to solve the system of 
first-order Dirac-Coulomb equation with the rough boundary condition or the exact boundary con- 
dition. We have not understood why many papers treating of the relativistic quantum mechanics 
seek the corresponding second-order Dirac equation for writing the so-called simplified solutions. 
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It should be pointed out that the solutions of the original system of first-order Dirac-Coulomb 
equation are not simplified, and one cannot obtain any simplified solution by constructing the 
corresponding second-order Driac-Coulomb equation, unless introducing some incorrect equations. 
We can use some mathematical theorems of the optimum differential equationsQQ to discuss 
such kind of problems [3], and find many papers for constructing second-order Dirac-Coulomb 
equations are incorrect in mathematical and physical signification. It must be ingeminated that 
the boundary condition and the uniqueness of solution are very important in solving wave equa- 
tion. Some classical {igI] [itI [l^ and modern[l^[2^[2l| quantum mechanics textbooks have treated 
these problems Not anyone should make any mathematical mistakes to spell backward the 



formula of the energy eigenvalues. The omnifarious mathematical mistakes concealed in the men- 
tioned paper are cross-sectional, implying there are too much similar problems in many published 
papers [2c 
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